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ABSTRACT 

We have a system of N sendox machines (service centers) 
which can transmit or receive (but not simultaneously) photo- 
copies to and from each other. The stochastic process 
(Xj^j(t),t > 0}, where Xj^j(t) is the number of messages under- 
going transmission in the system at epoch t, is selected to 
be the measure of system congestion. Under the assumptions 
of exponential inter-arrival and transmission times and that 
messages not transmitted immediately are lost, Xj^(t) is a 
birth and death process. We call the model under these as- 
sumptions the "loss” model. It is shown that Xj^(t) can be 
approximated by a diffusion process. The steady state re- 
sults of the diffusion model can be calculated with ease and 
are in close agreement with those of the birth and death 
model and the simulation "loss" model. A more realistic "no 
loss" model is also developed where messages that are not 
transmitted immediately are allowed to queue up. It is shown, 
however, that this "no loss" model can be quite accurately 
approximated by an appropriate "loss" model. 
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NOTATIONS 



N total number of machines (service centers) in the sys- 

tem. (Assumed to be an even number.) 

M=N/2 maximum number of messages that can be transmitted by 
the system. 

X arrival rate of messages at each machine. 

effective arrival rate of messages at each machine in 
the "loss model. 

M departure rate of a message undergoing transmission. 

Xj^(t) number of messages undergoing transmission at epoch t. 

Sw(t) stochastic element associated with X„(t). 

M M 

y(t) deterministic component of Xj^(t). 

(J)j^j(t) characteristic function of Xj^(t). 

i^mC^) characteristic function of Sj^(t). 

Xg steady state mean of Xj^^(t). 

steady state variance of Xj^(t). 

Pj ' s steady state probabilities of Xj^(t). 
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I. INTRODUCTION 



The main purpose of this thesis is to construct and com- 
pare various models which can be used to analyze the degree 
of congestion in a system of N sendox facsimile machines. 
Sendox machines can both transmit and receive (but not simul- 
taneously) photocopies using standard telephone lines. The 
original problem arose out of a proposal to install a number 
of these machines in central London Mintech buildings for 
transmission of thin "immediate" messages between buildings. 
This problem was first studied by Coleman [Ref. 1] who in- 
vestigated the expected system congestion, expected cost per 
message and a possible priority system. 

One of the models in his study considered each machine 
individually as a single server system. It was assumed that 
only one machine was located in each building. Arrivals at 
each of the N machines occur from two sources: 1) requests 
to transmit messages to other machines and 2) incoming mes- 
sages from the other N-1 machines. These two arrival streams 
were merged into one and the standard queuing formula for a 
single server queue was applied to obtain the expected wait- 
ing time of messages at an individual machine. However, a 
subsequent simulation study showed wide discrepancy between 
the theoretical and simulated expected waiting times. This 
occurred because the single server model did not take into 
account the "interference' between machines. The model as- 
sumed that when a message is transmitted from one machine, 
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its destination (which could be any one of the other N-1 ma- 
chines) is always available to receive the message. In 
reality, its destination may be engaged in transmitting or 
receiving another message. The effect of such "interference" 
tends to increase the waiting time as compared to the result 
obtained from the single server queuing model. Any reason- 
able theoretical model, therefore, has to consider all N ma- 
chines simultaneously. It is very difficult, however, to 
build an . analytical model that can give a reasonable esti- 
mate of waiting time because it has to consider N queues all 
at once. 
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II. DESCRIPTION OF THE "LOSS" MODEL 



The problem becomes much more tractable when we only con- 
sider modelling the system size i.e., the number of messages 
undergoing transmission. This number varies from 0 to M 
where M=[N/2], the largest integer less than or equal to N/2. 
Let Xjj(t) be defined as the number of messages undergoing 
transmission at epoch t. The stochastic process {Xj^j(t), 
t ^ 0} will be our measure of system congestion. 

A number of assumptions are made to facilitate construc- 
tion of theoretical models for the process Xj^(t). We assume 
a Poisson arrival of messages at each machine. The arrival 
rate of messages requesting transmission at each machine is 
the same and is denoted by X. Moreover, we assume that each 
message undergoing transmission is finishing at an exponen- 
tial rate of y. Furthermore we define p to be and hence 
the departure rate (y) of a message is also equal to 2pX. 
Finally, we assume that if a message fails to be transmitted 
immediately upon arrival (due to either the receiving or 
sending machine being busy), then the message is lost to the 
system. For this reason, we call the model based on the 
above assumptions the "loss" model. 

Two theoretical and two simulation models for the proc- 
ess will be described and discussed in the following 

sections. We are mainly interested in steady state behavior 
of the process Xj^(t). In the first theoretical model, Xj^j(t) 
is shown to be a birth and death process (this was done in 
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[Ref. 1]) whose steady state probabilities can be derived 
and- calculated numerically. Since Iglehart [Ref. 3] had 
proved the convergence of certain birth and death processes 
to diffusion processes, we were led to consider a diffusion 
approximation for the process Xj^(t), when M is large. This 
was the approach taken in the second theoretical model where 
it was shown that Xj^^(t) can be approximated by a non-station- 
ary Ornstein-Uhlenbeck process. Finally, a simulation model 
was built to validate the theoretical models. The steady 
state results of the theoretical and the simulation models 
compared very favorably with one another. 

A second simulation model was developed for the case in 
which messages not transmitted immediately, queue up until 
they are finally transmitted. We call this model the "no 
loss" model. We are unable to construct a theoretical model 
for this case. However we will show that this "no loss" 
model can be quite accurately approximated by the diffusion 
"loss" model. 
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III. BIRTH AND DEATH "LOSS" MODEL 



In an unpublished paper [Ref. 1], Dr. Rodney Coleman 

modelled the number of messages undergoing transmission, 

X,,(t), as a birth and death process. For the sake of sim- 
M 

plicity, we assume that N, the total number of machines in 
the system is even and hence M, the maximum number of mes- 
sages that can be carried by the system is given by N/2. 

When the process Xj^^(t) is in state j, the conditional proba- 
bility that Xj^(t) will visit state j+1 in the next interval 
of length At is X .At + o(At) and the conditional probability 
that will visit state j-1 is V^At + o(At), where X^ 

and vij are the infinitesimal birth and death rates respec- 
tively [Ref. 6], Messages can only enter the system at the 
(2M-2j) free machines at a rate of (2M-2j)X. When a message 
arrives at a free machine, its probability of finding an 
available receiving machine is > therefore, 



Since there are j messages, each finishing transmission at 
the rate of y, and thus 



X . 
3 



(2M-2J )(2M-2.i-l)X 
(2M-1) 



(3.1) 



y . = jy = 2jpX 



(3.2) 
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Let Pj denote the steady state probability of being in 

state j. Then Coleman [Ref. 1] showed that 



/ 



k! (2M-2k) ! 

M 

k=0 k! (2M-2k) ! 



k=0,l, . . .M 



where Y = ( 



1 

2M-1 




1 

2p(2M-l) 



However, since the caluclation for the P, 's are carried 
out on a digital computer, a more efficient method computa- 
tionally is to use the following recursive relationships 
[Ref. 6]. 

Pj^Pj^ = ^k— l^k— 1 k=l ,2,...M (3.3) 

M 



The computer program which calculates the steady state prob- 
abilities and also the mean and variance of steady 

state is listed in Appendix D. 
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IV. DIFFUSION "LOSS" MODEL 



We follow Gaver ' s [Ref. 4] approach by approximating the 

process Xj^j(t) with a deterministic component y(t) and a 

stochastic element or noise S,,(t). Thus X,,(t) can be written 

M M 

as : 



We show that as we let M go to infinity, Sj^(t) converges 
to a non-stat ionary Ornstein-Uhlenbeck process, S(t) [Ref. 

2] . Hence for large M, we can use the following formula to 
approximate Xj^(t): 



Because of our assumption of Poisson arrivals, a message 
will arrive at a machine with probability Adt + o(dt) in the 
time interval (t,t+dt). During this time interval, the 
probability of two or more messages arriving at the same ma- 
chine is o(dt). At epoch t the number of messages under- 
going transmission is Xj^(t) and hence the number of unoccupied 
machines is given by N-2Xj^(t) or 2M-2Xj^(t). In the time in- 
terval (t,t+dt), effective arrival of messages can occur 
only at any of the 2M-2Xj^j(t) unoccupied machines. Given that 
a message has arrived at an available machine, its probabil- 
ity of finding the receiving machine not in use is now 



Xj^j(t) = My(t) + /M Sj^(t) 



(4.1) 



Xj^j(t) = My(t) + /M S(t) 



(4.2) 



(2M-2Xj^j(t)-l)/(2M-l) . Therefore, 

Prob{X,,( t+dt )=X,,( t )+l 1 given X,,( t ) } 
M M ’ M 



(2M-2Xj^j(t))(2M-2Xj^j(t)-l)Adt 



(4.3) 



2M-1 



+ o(dt) 
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Similarly, under the assumption of exponential transmis- 
sion times, each message undergoing transmission is finish- 
ing at the rate of y or 2pX. At epoch t, messages will be 
finishing at a rate of 2Xj^^(t)pX. Therefore during the time 
interval (t,t+dt) the probability of a message completing 
its transmission is 2Xj^^(t)pXdt + o(dt) and the probability 
of two or more messages completing their transmission is 
o(dt). Therefore, 

Prob {Xj^(t+dt) = Xj^(t)-l| given Xj^(t)} 



= 2X^pXdt + o(dt) 



(4.4) 



Finally, 

Prob (Xj^j(t+dt) = X^^(t)| given X^(t)} 



M’ 



= 1 - 



(2M-2Xj^(t))(2M-2Xj^(t)-l) 

2U^ 



M' 

(4.5) 

Xdt - 2Xj^(t)pXdt + o(dt) 



In Appendix A, equations (4.1) through (4.5) are used to 
derive the partial differential equation 



dt 



-2X{2+p-2y(t ) }e-|^ 



2X{y2(t)-y(t) (2-p)+l}|-i(j 

^ (4.6) 



Here ij^(6,t) is the characteristic function of S(t), the lim- 
it of the noise process Sj^(t) as M ^ The deterministic 
process y(t) must satisfy the ordinary differential equation 



y’(t) - 2X{y2(t) - (2+p)y(t) - 1} = 0 (4.7) 

in order that the partial differential equation (4.6) be 
valid. The partial differential equation (4.6) is recognized 
as the transformed version of the forward differential 



17 



equation for a non-stat ionary Ornstein-Uhlenbeck process. 

It is known [Ref. 2] that the probability density function 
of the Ornstein-Uhlenbeck process S(t) must satisfy the 
Kolmogorov Forward Differential Equation 

If = {a(s,t)p} - 1^ {B(s,t)p} (4.8) 

where $(s,t) and a(s,t) denote the infinitesimal mean and 
variance respectively, conditioned on S(t)=s. 

3(s,t) = -2X{2+p-2y(t)}s (4.9) 

a(s,t) = 2X{y2(t)-(2-p)y(t)+l} (4.10) 



The first step toward determining the mean and variance 
of the process S(t) involves solving the first order differ- 
ential equation (4.7) to obtain y(t) explicitly, assuming 
the initial condition that y(o)=0. The solution is derived 
in Appendix A and is 



y(t) 



= v-w { 



1 + ke 



-4Xwt 



-} 



1 - ke 



-4Xwt 



(4.11) 



where v = 1 + ip 

= (1 + ip)2 - 1 

k = 

w+v 



Since the Ornstein-Uhlenbeck process is Gaussian, we may 
differentiate the Gaussian characteristic function 
exp [i6a(t )-i0^b(t ) ] with respect to 6 and t, substitute ip , 
■|^, back into equation (4.6) and equate the coefficients 
of i0 and 0^ in order to obtain two first order differential 
equations for a(t) and b(t) which denote the mean and 
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t 






variance of the process S(t) respectively. Assuming that 
S(0)=0, we can show that 



E[S(t)] = 0 



t > 0 



(4.12) 



and since 



Xj^(t) - My(t) + /M S(t) 



we obtain 



E[Xj^(t)] - My(t) t > 0 



(4.13) 



where y(t) is given by equation (4.11) 

The solution of the first order differential equation for 
the Var {S(t)> involves integrals that are quite intractable 
and we are only able to obtain an approximation for large t. 



yi " yp(2-p)+i 
l(2fp-2y^) " 



(4.14) 



where y^ = y(t) = 1 + ip - /(l+ip)^-l 



(4.15) 



Since Var {S(t)> - Var {X^(t)>, we obtain the following 
expressions for steady state mean and variance of Xj^(t), de- 
noted by X and respectively, 
s s 



X = E [ X. . ( t ) ] = My 
s ^ •'e 

ol = Var {Xj^j(t)} 



t large 



(4.16) 



« l 2T^2^) -l 



t large 



(4.17) 



Detailed derivation of equations (4.12) and (4.14) is given 
in Appendix C. When t is large, X,,(t) is normally distri- 
buted with moan and variance given by equations (4.16) and 
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! 

I 



(4 . 17) respectively . Plots of x /M and a^/M against P are 

s s 

given in Figures 1 and 2 respectively. It is seen in Appen- 
dix C that y(t) reaches its steady state value of y^ consid- 
erably- faster than Var {S(t)}, thus the steady state mean of 
X,,(t) is attained much faster than the steady state variance 
of 



The steady state probabilities for X,.(t) can be easily 

M 

obtained from the standard normal table after x and are 

s s 

calculated using equations (4.16) and (4.17). However, since 
we are using a continuous (normal) probability distribution 
to approximate a discrete distribution, we have to make the 



following "continuity corrections" 

2 



i 

O — oo 



/^o 



e dx 



i+4 

■ P. = ./i — — 



(X-Xg) 
— 



M 



00 

= / 



/2na, 



2o‘ 



(x-x^) 



dx 



dx 



i=l,2, . . .M-1 



(4.18) 



(4.19) 



(4.20) 



Another way to derive the parameters for the Ornstein- 
Uhlenbeck process S(t) is by defining the stochastic differ- 
ential equation for the process Xj^(t). 

dXj^j(t) = 3j^(x,t)dt + Z(t) /ttj^(x, t )dt (4.21) 

where Z(t) is a purely random Gaussian process with zero 

mean and unit variance and 3j^^(x,t)dt and Uj^(x,t)dt are the 

mean and variance of the increment dX,,(t) in a small time 

M 
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interval (t,t+dt). Recall that we have derived the follow- 
ing probabilities earlier 

Prob (dXjj(t) = l| given Xj^(t)} 



{ 2M-2Xjj( t ) } { 2M-2Xj^( t )-l) 

2M-1 



Xdt + o(dt) 



Prob {dXj^j(t) = -l| given Xj^(t)) 
= 2Xjj(t)pAdt + o(dt) 



Prob {dXj^(t) = 0| given Xj^(t)} 

{ 2M- 2X ( t ) } { 2M- 2X ( t ) - 1 } 

= 1 ^ Xdt - 2Xj^(t)pXdt + o(dt) 



Therefore it is reasonable to define 3j^(x,t)dt and aj^(x,t)dt 
as follows: 



3j^j(x,t)dt = E[dXj^(t)] 

{2M-2Xj^(t)}{2M-2Xj^(t)-l} 

2M-1 



(4.22) 

Xdt-2Xj^(t)pXdt 



aj^(x, t )dt 



Var [dXj^(t)] 

{2M-2Xj^(t )}{2M-2Xj^j(t)-l) 
2M-1 



(4.23) 

Xdt+2X„(t)pXdt 

M 



Substituting Sj^(t) as given by equation (4.1) in place of 
XM(t), we derive a stochastic differential equation for 
Sj^j(t) in Appendix B. Then by letting M and requiring 

the first order differential equation (4.7) for the deter- 
ministic process y(t) to hold, we obtain the stochastic 
differential equation 
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dS(t) = 2X{2y(t)-2-p}S(t)dt 



+ Z(t) /2X{y^ (t)-y(t)(2-p)+l}dt (4.24) 

This equation is the stochastic differential equation for an 
Ornstein-Uhlenbeck process which has the same infinitesimal 
mean and variance as those given in equations (4.9) and (4.10). 
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V. SIMULATION "LOSS” MODEL 



A. SIMULATION PROGRAM I 

Two methods of simulations are presented here for the 
"loss" model. Simulation Program I makes use of the obvious 
approach of duplicating the operation of the actual system. 
Messages arrive at each of the N machines at a Poisson rate 
X. These messages are either lost or undergo transmissions 
which are being completed at an individual rate of y. A 
brief description of the program logic is given in the next 
paragraph. 

Simulation Program I is written in FORTRAN. It operates 
simply by searching for the next event to occur. In this 
model an event is either the arrival of a message at a ma- 
chine or the completion of a message undergoing transmission. 
Each event is assigned a clock time which denotes the time 
of occurrence of that event. The next event is that whose 
clock time is the minimum of the clock times of all the fu- 
ture events. Having found this event, the program performs 
a series of tests to determine whether the conditions in the 
system are such that the potential event can take place. If 
so, the program moves to execute this event. For example 
the event may be the arrival of a message at machine i, then 
the necessary test is to check whether the sending machine 
i and a randomly selected receiving machine j , j^i are both 
available. If so, the number of messages undergoing trans- 
mission increases by one. The basic loop is repeated as 
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many times as necessary to complete the entire simulation. 
The flow chart and listing of this program are given in Ap- 
pendix D. 



B. SIMULATION PROGRAM II 

Another way to develop a simulation program for the "loss" 



model is to make use of the fact that X,,(t) is a birth and 

death process. We know the process Xj^^(t) stays in state i 

for an amount of time that is exponential with rate 

and given that it makes a transition from state i, it goes 

y i 

to either state i-1 or i+1 with probabilities -r — r — and 
— T — respectively where X. and y. are as defined in section 

A i"'‘y i 1 1 

III. 



Instead of monitoring the flow of messages through the 
system of N machines,' Program II works by keeping track of 
the length of time that the process Xjj(t) spends in each 
state and when Xj^(t) changes state. This program is also 
written in FORTRAN. It starts off with Xj^(O) in a given 
initial state i. An exponential variate with mean ^ ^ ■ 
is generated and this represents the length of time that 
Xj^(t) spends in state i before jumping to state i-1 or i+1. 
In order to decide whether Xj^(t) next visits i-1 or i+1, a 
uniform random variate p, 0 _< p _< 1, is generated. If p is 
less than -r — - — , then X (t) visits state i+1 next and if p 



Xi+yi 



M 



exceeds 



-r — - — X,,(t) goes into state i-1. Once X„(t) switches 
Xj[+yj M M 



to a new state j, we repeat the process of determining the 



length of Xj^j(t)'s stay in state j and the next state to be 
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visited. If Xj^(t) is in state 0, it will next visit state 
1 with probability one. Similarly if Xj^j(t) is in state M, 
it will next visit state M-1 with probability one. The out- 
put of this program will be discussed together with that of 
Program I in the next section. The flow chart and listing 
of Program II are included in Appendix D. 

C. DISCUSSION OF OUTPUT OF PROGRAMS I AND II 

Both simulation programs have five input parameters, 
namely the total number of machines (N), the duration of each 
simulation run (TSTOP), the arrival rate of messages at each 
machine (X), the departure rate of transmitted messages (u) 
and the initial state i. Since we are mainly interested in 
steady state results, the primary output of both simulation 
programs is the probability distribution of the process 
XM(t) when in steady state. 

One way to obtain these steady state probabilities is to 
let the simulation program run for t' minutes until we are 
reasonably sure that steady state has been achieved, and 
then record that value of X,,(t'). Since each simulation run 
produces only one realization of Xj^j(t'), this method of es- 
timating the P.'s proves to be very inefficient because of 
the large number of simulation runs required to produce a 
consistent estimate. 

A better method to estimate the P.'s is to observe the 

J 

time average probabilities of X,,(t) when it is already in 

M 

steady state. We choose a starting state i which is very 

close to the steady state mean My so that X,,(t) will reach 

e M 
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steady state shortly. The simulation model will have to run 
for a considerable amount of time in order to gather enough 
sample values. However, this results in a program that re- 
quires too much core storage. To overcome this problem, we 
made 200 simulation runs of TSTOP minutes instead of a single 
run of 200 x TSTOP minutes. 

During each execution of the simulation program, 200 
simulation runs of duration TSTOP minutes each are made. 

Let fjj^ denote the frequency that the process Xj^(t) visits 
state i. At the beginning of the program, all the are 

set equal to zero. During each simulation run, the number 
of messages undergoing transmission is observed at discrete 
one minute intervals. If this number is k, then fj^ is incre- 
mented by one and all other f.'s, j^^k, remain unchanged. 

Thus, during the entire simulation program, Xj^(t) is observed 
at 200 X TSTOP points, and the estimated steady state proba- 
bilities are given by: 

- = ^k 

Pk 200XTSTOP 

Initially, both programs are tested with the same input 
parameters to ensure that they are yielding results that are 
consistent with each other. From then on, only Program II 
is used to generate the various output that are used to val- 
idate the theoretical models, because time-wise, it is much 
more efficient than Program I. For example, with N=50, 

A=2.5 and y=10.0, the total execution time of Program I is 
about 18 times that of Program II. 
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VI. SIMULATION "NO LOSS" MODEL 



In the previous sections, we have developed analytical 
and simulation models for a system of N machines in which 
messages that are not transmitted immediately are lost. A 
more realistic assumption is to allow these messages that 
are not transmitted immediately to join a queue at the re- 
spective machines. We call a model based on this less strin- 
gent assumption a "no loss" model. This "no loss" model is 
not easily amenable to theoretical analysis, but it can be 
easily simulated by a program written in GPSS. The assump- 
tions of exponential distribution for the inter-arrival and 
inter-departure times of messages could now be relaxed be- 
cause a simulation program can easily handle any other prob- 
ability distribution. However, the previous assumption of 
exponential inter-arrival and transmission times is retained 
for comparison purposes. 

In the GPSS program, the Sendox machines are treated as 
facilities and the messages that arrive at the various ma- 
chines are treated as transactions. The program consists 
mainly of N closely similar segments. Each segment simulates 
the activity of a facility. The program instructions for 
each segment are almost identical except for the fact that 
each facility has a different facility number i, i=l,2,...N. 

Arrivals of messages are simulated through the generation 
of random arrivals of transactions at an exponential rate X 
at each facility. A transaction created at facility i is 
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randomly assigned to a destination facility j, inde- 

pendently of other transactions. Then facilities i and j 
are tested to see whether they are available simultaneously. 

If so, the transaction proceeds to seize both facilities for 
a length of time that is exponentially distributed with rate 
M. If not, the transaction will have to wait until both 
facilities i and j are simultaneously free. If another trans- 
action is created while a previous transaction is either 
waiting for or undergoing service, the new transaction must 
queue up. Therefore we assume a first come first serve 
queue discipline. This procedure is followed at all N facil- 
ities. 

GPSS is preferred to FORTRAN in this model because the 
logic and structure of a GPSS program are more simple. More- 
over queue statistics and steady state probabilities of 
XM(t) can be gathered with little extra effort. The only 
shortcoming is that this program cannot adapt easily to 
changes in the number of machines. Thus if we want to double 
the number of machines the card deck will almost double its 
size. This GPSS program was only written for N equal to 20, 
whereas N can go as high as 200 in Programs I and II. A 
listing of this program is included in Appendix D. 
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VII. COMPARISON OF RESULTS 



According to the diffusion "loss” model, when Xj^(t) is 

in steady state Xj^^(t) is approximately normally distributed 

with mean x and variance given by equations (4.16) and 
s s 

(4.17). Moreover, using x and and the continuity correc- 

s s 

tion, given by equations (4.18), (4.19) and (4.20), the 
steady state probabilities can be obtained from the standard 
normal table. 

We are mainly interested in determining the accuracy of 
the diffusion "loss" model for various values of N and p in 
comparison with the simulation models. For ease of refer- 
ence equations (4.16) and (4.17) are repeated here. 



= My 



N 

2 ^e 



= M 



y|-ye(2-P)+i 

2(2+p-2y^) 



(7.1) 

(7.2) 



where y^ = 1 + ip - /(l+4p)^-l 

Equations (7.1) and (7.2) implies that x and depends 

s s 

only on the ratio of X and y , when M is kept constant. To 

check this implication. Simulation Program II was run with 

five sets of X and y where the ratio is kept equal to 2. 

The steady state means and variances obtained from these 

runs turned out to bo extremely close to x and as calcu- 

s s 

lated from equations (7.1) and (7.2) and these results are 
given in Table I. 
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Table I 



Steady State Mean and Variance of with 



p constant, N 



100, p = fx " 



Arrival rate 
of messages 

X 

per hour) 


Departure rate 
of messages 

y 

per hour) 


Simulation 
"loss" model 


Diffusion 
"loss" model 


X 

s 


s 


X 

s 


s 


1.25 


5.0 


13.389 


7.712 


13.389 


7.723 


2.5 


10.0 


13.406 


7.753 


13. 389 


7. 723 


5.0 


20. 0 


13.418 


7.787 


13.389 


7.723 


10.0 


40.0 


13.410 


7.760 


13.389 


7.723 


20.0 


80.0 


13.408 


7.760 


13.389 


7.723 



(Each simulation runs for 2400 hours) 
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Next, keeping M=10, or N=20, six more runs were made with 

the values of p ranging from 0.5 to 10. Again x and 

s s 

calculated from equations (7.1) and (7.2) were in close agree- 
ment with the simulation results. These results are summa- 
rized in Table II. 

Equations (7.1) and (7.2) also imply that x and are 

s s 

linear functions of N keeping p constant at 2. Six simula- 
tion runs were made for values of N ranging from 10 to 200. 

The steady state probabilities for the birth and death "loss" 
model were also calculated for the same values of p and N. 
Again, the results indicated that all three models agree very 
well as far as steady state mean, variance and cumulative 
probabilities are concerned. The steady state means and 
variances for the three models are given in Tables III and 
IV respectively. In Table V we compared the cumulative 
probabilities of the three models for the case where N=20 
and p=2. These same values when plotted on normal probabil- 
ity graph paper (Figure 3), gave clear visual evidence of 
their close agreement. 
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Table II 



Steady State Mean and Variance of with 

N constant, N = 40, X = 2.5. 



Departure rate 
of messages 

y 

(# per hour) 




Simulation 
"loss" model 


Diffusion 
"loss" model 


X 

s 


a" 

s 


X 

s 


s 


2.5 


0.5 


9.976 


3.647 


10.000 


3.333 


5.0 


1.0 


7.631 


3.544 


7.600 


3.413 


7.5 


1.5 


6.263 


3.357 


6.277 


3.299 


10.0 


2.0 


5.352 


3.117 


5.359 


3.094 


30.0 


6.0 


2.537 


1.969 


2.540 


1.967 


50.0 


10.0 


1.676 


1.421 


1.678 


1.415 



(Each simulation runs for 2400 hours) 
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Table III 



Steady State Mean of Xj^(t) with X and y constant, 
X=2.5, y=10.0 and p = = 2.0. 



Number of 
machines 
N 


Steady State Mean, 


Simulation 
"loss" model 


Diffusion 
"loss” model 


Birth and Death 
"loss” model 


20 


2.703 


2.680 


2.697 


40 


5.373 


5.359 


5.376 


60 


8.049 


8.038 


8.055 


80 


10.722 


10.718 


10.735 


100 


13.397 


13.398 


13.414 


200 


26.803 


26.795 


26.811 



(Each simulation runs for 2400 hours) 
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Table IV 



Steady State Variance of Xj^^(t) with X and y 
constant, X = 2.5, y = 10.0 and p = = 2.0. 



Number of 
machines 
N 


Steady State Variance 

s 


Simulation 
"loss” model 


Diffusion 
"loss" model 


Birth and Death 
"loss" model 


20 


1.547 


1.547 


1.559 


40 


3.052 


3.094 


3.106 


60 


4.617 


4.641 


4.653 


80 


6.142 


6.188 


6.199 


100 


7.723 ■ 


7.735 


7.746 


200 


15.440 


15.470 


15.484 



(Each simulation runs for 2400 hours) 
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Table V 



Steady State Cumulative Probabilities, N=20, 
M=N/2=10, A = 2.5, u = 10.0, P = ^ =2.0. 





Steady State Cumulative Probabilities P^ 


State j 


Simulation 
"loss” model 


Diffusion 
"loss” model 


Birth and Death 
"loss" model 


0 


0.0275 


0.0401 


0.0280 


1 


0.1646 


0.1716 


0.1681 


2 


0.4459 


0.4427 


0.4501 


3 


0. 7485 


0.7453 


0.7470 


4 


0.9243 


0.9284 


0.9247 


5 


0.9865 


0. 9887 


0.9864 


6 


0.9985 


0.9989 


0.9986 


7 


0 . 9999 


1.0000 


0.9999 


8 


1.0000 


1.0000 


1.0000 


9 


1.0000 


1.0000 


1.0000 


10 


1.0000 


1.0000 


1.0000 



(Each simulation runs for 2400 hours) 
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Figure 3 
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CUMULATIVE PROBABILITY (xlOO) 



VIII. APPROXIMATION OF "NO LOSS" MODEL BY "LOSS" MODEL 



Since the diffusion "loss" model gives relatively simple 
formulas for x , and P.'s, it appeared desirable to re- 
late the results of the "loss" model to the "no loss" model. 
It may be noted that in the "loss" model, only a fraction of 
the arriving messages are transmitted, whereas in the GPSS 
"no loss" model all arriving messages are eventually trans- 
mitted. 

First, we need to determine the fraction of messages 
transmitted in order to obtain the effective rate of message 
arrivals at a "loss" model. In the "loss" model, the prob- 
ability that a message arrives at a busy machine is approxi- 
mately X divided by M or y , assuming steady state. Since 

S G 

a message can be transmitted only if both the transmitting 
and receiving machines are available, the probability of a 
message being transmitted immediately is approximately 
(1-y^)^ assuming independence between availability of ma- 
chines and that N is large. Hence is the effective 

arrival rate of messages at a machine in the "loss" model is 
given by 

Xeff = Ml-Yg)' (8.1) 

This leads to the idea that a "no loss" model with input 
parameters N, X and y can be approximated by a "loss" model 
with input parameters N, X' and y where X' is such that 
X'(l-y ' )^ = X. Using equation (4.15) for y ' with p' = 

G G ^ A 
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the equation X'(l-y^' )^ = X can be rewritten as 




( 8 . 2 ) 



After more rearrangement of the various terms (details are 
given in Appendix E), we obtain an equation expressing X' 
explicitly as a function of y and X: 



The above relationship between "loss" and "no loss" mod- 
els was tested using the simulation models with five differ- 
ent values of y and X' as input parameters, N and X being 
kept constant at 20 and 2.5 respectively. The steady state 
means and variances are summarized in Table VI. For the 
case of N=20, y=10.0, X=2.5 and X'=10.0, the steady state 
cumulative probabilities are also calculated using the dif- 
fusion "loss" model, and the results from all three models 
are tabulated in Table VII and plotted in Figure 4. All the 
results indicate that the approximation of the "no loss" 
model by the "loss" model is quite accurate. 



X 



X ' 



( 1 - 



(8.3) 
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Table VI 



Comparison of "no loss" model with its 
"equivalent" "loss" model, N=20, A=2.5. 



Simulation "no loss" 


model 


Simulation 


"loss model 


X 


y 


X 

s 


a" 

s 


X’ 


^eff 


X 

s 


a" 

s 


2.50 


7.50 


6.677 


1.362 


22.500 


2.50 


6.694 


1.343 


2.50 


8.57 


5.758 


1.503 


14.407 


2.50 


5.862 


1.556 


2.50 


10.00 


4.972 


1.640 


10.000 


2.50 


5.031 


1.659 


2.50 


12.00 


4.105 


1.732 


7.347 


2.50 


4.192 


1.726 


2.50 


15. 00 


3.319 


1.651 


5.625 


2.50 


3.357 


1.680 



(Each simulation runs for 2400 hours) 
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Table VII 



Steady State Cumulative Probabilities Calculated 
from Three Different Models, N=20, y=10.0. 



State 

3 


Steady State 


Cumulative Probabilities P. 

3 


Simulation 
"no loss" model 

X=2.5 


Simulation 
"loss" model 

X'=10. 0 


Diffusion 
"loss" model 

X'=10.0 


0 


0. 0000 


0.0001 


0.0002 


1 


0.0010 


0.0028 


0.0034 


2 


0. 0240 


0.0245 


0.0266 


3 


0. 1210 


0.1159 


0.1146 


4 


0. 3480 


0.3347 


0.3492 


5 


0.659 


0.6380 


0.6507 


6 


0.8900 


0.8780 


0.8874 


7 


0. 9810 


0.9790 


0.9734 


8 


0.9980 


0.9985 


0.9966 


9 


1.0000 


0.9999 


0.9998 


10 


1.0000 


1.0000 


1.0000 



(Each simulation runs for 2400 hours) 
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IX. CONCLUSION 



It was found that the steady state results of the diffu- 
sion and the birth and death models were in very close agree- 
ment with the simulation model, through comparison of the 
steady state means, variances and cumulative probabilities. 
The usefulness of the diffusion model, as compared to the 
birth and death model, lies in the ease with which its steady 
state mean and variance and its discretized probabilities 
can be computed without resorting to a digital computer. 

The simple form of equations (7.1) and (7.2) exhibits 
the fact that the steady state mean and variance vary line- 
arly with N and depend only on p, and not on the individual 
values of X and y. These results cannot be inferred from 
the steady state formulas of the birth and death model or 
the simulation model. 

Lastly, we developed a simulation model for the much 
more complicated case when the messages are allowed to wait 
for transmission. We showed that a "no loss" model with in- 
put parameter N, X and y can be quite accurately approxi- 
mated by a "loss" model with input parameter N, X', and y 
where X' is calculated through equation (8.3) from the other 
parameters. 
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APPENDIX A 



Derivation of the Characteristic Function of S(t) 



Define <f>jj(6,t) and the characteristic func- 

tions of Xj^^(t) and Sj^(t) respectively: 



= E[e 



iexM(t) 



ies (t) 

iPj^(0,t) = E[e ] 



(A.l) 
(A. 2) 



Since 



Xj^j(t) = My(t) + M Sj^(t) 



we obtain 

I^j^(6,t) = t) (A. 3) 

Given that Xj^(t) = j, equations (4.3), (4.4) and (4.5) give 
the respective probabilities that Xj^^(t+dt) be equal to j+1, 
j-1 or j . Hence we can write down the following expression 
for 4>j^j( 6 > t+dt ) : 

(})jj(e, t+dt) = 

= E[e^®^^M^''^^'^^^{4M2-2M+2Xj^(t)(l-4M)+4Xj^j(t)}||^ 

+ e^®^M^^^{l-2Xj^j(t)pAdt-(4M2-2M+2Xj^j(t)(l-4M) 

e^®^^M^^^"^^{2Xj^(t)pAdt) + o(dt)] 

(A. 4) 
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Then after rearranging the terms and taking the limit as 
t we obtain 

^ 4>j^(0 , t+dt)-(})j^(6 , t) 

at dt-^0 dt 



= E[ e^®^^^^^^e^®{4M^-2M+2Xj^j(t)(l-4M)+4X2 (t)>2^3j 
+ e^®^M(^){_2pXXjj(t)-(4M^--2M+2Xj^(t)(l-4M)+4Xj|(t)) 

+ e^®^M(t)g-ie^2pXXjj(t)}] (A. 5) 

Further rearranging of terms results in 



at Ke ■*•>>211-1 



•f 2 XXj^^(t)e^®^M<^^){|jj^(e^®-l)+p(e“^®-l)} 

+ 2XMe^®^^^^'^\e^®-l)] (A. 6) 



However, since 



94 >ji(e,t) 

ae 



= E[iXjj(t)e^®^M^‘‘^^] 



aH^Ce.t) 



= E[-X^(t)e^®^M(‘t)] 



the above can now be written as 



9<}>^^(e,t) ^ ie 4X ^ ie .. 

at aF^^ ® ^ 2 M-i ^ 2 XMc})j^j(e,t)(e i) 



+ 2Xi 



_94>fj(6,t) A^_i ie- ^ , -ie 



ae 



{ 



2M-1 



(e" -l)-p(e " -1)} 



(A. 7) 
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Rewriting equation (A. 3) we obtain the following expressions 



for 4>,,(6/»^,t) and its first and second partial derivatives 

M 



<t.^^(6//M,t) = e 



_ i0/My(t) 






(A. 8) 



9 ^ = e U 3 ^^ + i6/My ' (t)ipj^j(e , t)} (A. 



9) 



^ ^ie/MyCt)^^®''^) 



90 



gg + i/My(t)ipj^(0,t)} (A. 10) 



92(|,j^(e//M,t) 

9lP 






90 



-My^ (t)if)j^(6 , t)> 



(A. 11) 



If we rewrite equation (A. 7) with 0//M in place of 0, we 
obtain 

94>j^(0//M,t) 

9t 

+ 

+ 2XM4>j^j( 0//M, t)(e^®/‘^-l) (A. 12) 



924>jj(e//M.t) 



2Av^i 



9cf)j^(0//M,t) 



,, 10//M,4AM 

>2MIT 



96 



r4M-l, iO//M -10//M 

^2ira^® -1)-P(e ' -1)} 



Using equations (A. 8) through (A. 11) to replace all the terms 
involving 4> by and factoring out equation 

(A. 12) becomes 
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3*„(6,t) 
3 1 



+ i0/My' ,t) 



= * 2i/My(t 



) My2 (t )({; ,(0 , t)} 

30 ” 



'M 9 i|j C 0 , t ) 



. 10//M -i0//M^Tr 

+ 2 X/Mi{ 2 j;pj(e ' -l)+p(l-e ' )}{ 



+i/My(t),|^jj(0,t)} + 2XM(e^®/’^“-l),|;jj(0,t) 



(A. 13) 



Substituting e with the 



ir power series 



^ _ gl0/»/M ^ _10 + |i + 0(-^) 

M3/2 



/M 



1 _ e = li + ii + o(-ir-) 

/M 



and after more rearranging of terms in (A. 13) we may write 



3i|; 



M 



3t 



+ i0/Miijj^{y' (t )-2X(7yi^y^ (t)-( 



^M-l-' 



4M-1 

2M-1 



+P)y(t)-1)} 



9i|^, ^ 



4XM . -i0 . 0" . r_4M_^ i0" . 

2M>-30^ 2/m 



^4M-1,,„, i0\, , r 2M 2,^^ 

- <2Mrr)«" 

- <lri -P)y(t)+i)2x(|i)t„+o(-^) 

/M 



(A. 14) 



where the arguments of (j^j^j(0,t) have been dropped for ease of 
exposition. The deterministic process y(t) must satisfy 
the following first order differential equation, as we let 
M ->■ oo 



y'(t) -2X {y2(t) - (2+p)y(t) -1} = 0 



(A. 15) 
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in order that equation (A. 14) simplify to 



= -6-||{2X(2+p-2y(t)}- |-i|; { 2X(y ^ ( t )- ( 2-p )y ( t )-l ) } 

(A. 16) 

This is the transformed version of the forward differential 
equation for a non-stationary O.U. process. 

We have to solve the first order differential equation 
(A. 15) in order to obtain y(t) as explicit function of t. 

For the sake of convenience, y'(t) and y(t) are written as 
^ and y respectively. Then equation (A. 15) is 

^ = 2X{y2-y(2+p)+l} (A. 17) 

Now let 

V = 1 + 4p 

and 

w^ = (1+ip)^ - 1- 

Then equation (A. 17) can be written as 

^ = 2X{(y-v)2-w^} 

which in turn may be written as 

( — - — - — )dy = 4Xwdt 
y-v-w y-v+w'^ 

Integrating both sides of the equation we find 



In 



y-v-w 

y-v+w 



4Xv/t + Cl 



or 



y-v+w 

y-v-w 



C 2 e 



-4Xwt 
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Assume that the system starts off empty, hence y=0 when t=0, 
then C 2 is evaluated as 

_ w-v 
w+v 



Now let 



Then 



w-v -4Xwt y-v+w 

^ rr: ■ - . 0 rr: ■ — - 

w+v y-v-w 



y = V - w 



Rewriting y as a function of t again 



y(t)=v-w{ 



- , w-v -4lwt 

1 + — ; — e 

w+v • 

w-v -4Xwf 

^ 



Further, 






lira y(t) 

t->oo 



= v-w 



(A. 18) 



(A. 19) 
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APPENDIX B 



Derivation of Stochastic Differential Equation 



From equations (4.21), (4,22) and (4.23) 



dXjy(t) = e^(x,t)dt + Z(t)/aj^(x, t)dt 



{2M-2X„(t)}{2M-2X (t)-l) 

B^l(x,t) = M ^ -2X^(t)pX 



{ 2M-2X ( t ) } { 2M-2X„( t )-l} 

ajj(x.t) = M ^ +2Xj^(t)pX 



(B.l) 

(B.2) 

(B.3) 



where Z(t) is a purely random Gaussian process with zero 
mean and unit variance. Using equation (4.1) we can express 
X,,(t) in terms of S,,(t) 

Xj^(t) = /MSj^(t) + My(t) 
therefore 

dXj^(t) = /MdSj^^(t) + My'(t)dt (B.4) 



Using equations (B.2), (B.3) and (B.4) to substitute for 

3(x,t), a(x,t) and dX,,(t) in equation (B.l) we obtain 

M 



AidSy(t) + My' (t)dt 

{2M-2/MS ,(t)-2My(t)}{2M-2v^S (t)-2My(t)-l) 

= I 2M3I 

- 2pX{*^Sj^j(t)+My(t)}dt] + Z(t)I2pX{/MSy(t)+My(t)}dt 



{2M-2/MS (t)-2My(t)}{2M-2»/MS (t)-2My(t)-l) , 

2MTI 
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Carrying out all the necessary multiplication, rearranging 
some of the terras and (since we are eventually going to let 
M go to infinity) approximating 2M-1 by 2M, we obtain the 
following expression that is exact in the limit. 



/MdS,,(t) + My'(t)dt 

M 

_ 1 

== ldt[ 2M-4My(t )+2My^ (t )+Sj^j(t ){4/M(y(t )-!)+ 



+2Sj^j2(-t-)_2p{/MSj^(t )+My(t)}] + Z(t)/Xdt [ 2M-4My(t )+2My2 (t) 



+Sj^(t){4/M(y(t)-l)+-;^}+2Sj^2(t)+2p{/MSj^(t)+My(t)}] 



Dividing by /M throughout the whole equation and after more 
rearranging of terms, we have 



dSj^(t) + /My’ (t)dt 

2/MXdt{y2(t)-(2+p)y(t)+l}+2XSj^(t)dt{2y(t)-2-p+ 

21 

+ + Z(t)[2Xdt{y2(t)-(2-p)y(t)+l} 

2XSjj(t)dt 



/M 



■(2y(t)-2+p+ ij) + 



(B.5) 



Now we let M go to infinity and assume that y'(t) = 2X{y^(t) 

-(2+p )y(t )+l} . Then with S(t) = lim S (t) we have' 

M^-oo ^ 

dS(t) = 2X{2y(t)-2-p}S(t)dt+Z(t)[2X{y^(t)-(2-p)y(t)+l}dt]^ 



Again we have shown that the process S(t) is an Ornstein- 
Uhlenbeck process with the following infinitesimal mean and 
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conditioned on S(t)=s. 



variance, 

3(s,t) = -3(t)s = -2X{2y(t )-2-p}s 

a(s,t) = a(t) = 2X{y^(t)-y(t)(2-p)+l} 
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APPENDIX C 



Derivation of Steady State Mean and Variance of Xj^(t) 



From equation (4.6) we have 

-|^ = -3(t)6-|^ - a(t)(-|— )\|^ (C.l) 

where 

3(t) = 2X{2+p-2y(t)} 

a(t) = 2X{y^ (t)-y(t)(2-p)+l) 

Let the mean and variance of S(t) be denoted by y(t) and 
x(t), and since S(t) is normally distributed, we can write 
its characteristic function \|;(6,t) in terms of y(t) and 
T(t ) . 

iK9,t) = (C.2) 

Differentiating i|;(6,t) with respect to 6 and t then dropping 
the arguments of i|^(6,t), we obtain 

ll = {iey ' (t)-ie2T* (C.3) 

= {iy(t)-0T(t))e^^®^^''^^"^®^'^^’*^^^ (C.4) 

O D 

We substitute these equations for tj;, into equation 

(C.l). Cancelling out the factor f('t)} 

tain 

ioy’ (t)-^e^T' (t) = -6(t)e{iy(t)-0T(t)}-ie^a(t) (c.s) 
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By equating coefficients of i6 and 6^ on both sides, we ob- 
tain the following two first order differential equations: 

p'(t) = -g(t)y(t) (C.6) 

T'(t) = -2B(t)T(t) + a(t) (C.7) 

Equation (C.6) can be easily solved by separation of vari- 
ables [Ref. 7], yielding 

t 

- / 3 ( X ) dx 

p(t) = y(0)e “ (C.8) 

It is reasonable to assume that the noise process S(t) starts 
off at time zero with zero mean and variance, hence y(0)=0 
and t(0)=0. 

Substituting the initial condition y(0)=0 into equation 
(C.8), we obtain 

y(t) = 0 for all t 0 (C.9) 

Equation (C.7) can be solved [Ref. 7] by the use of an inte- 
grating factor, giving us the following expression for x(t), 
with t ( 0)=0 

t z 

-2 f 3(x)dx t 2 / 6(x)dx 

x(t) = e ® {/ a(z)e ® dz} (C.IO) 

0 

Note that B(t) and a(t) are functions of y(t) where y(t) is 
a complex function of t, given by equation (A. 18): 



y(t) 



v 



w 



1 + 
{ — 
1 - 



w-v 

w+v 

w-v 

w+v 



-4Xwt 

e 

-4Xwt ^ 
e 



(C.ll) 
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Therefore the two integrals on the right hand side of equa- 
tion (C.IO) are quite intractable, unless t is very large, 
in which case we are able to derive an approximation for 
T(t). 

For the typical values of A and p that we are using, the 
“4 X wt 

term^ e in equation (C.ll) vanishes very rapidly with t 

approaching infinity. For example, with A=2.5 messages/hour, 
p=2 and t=1.0 hour, = 3.0 x 10“ ® . 

So for time t greater than one hour, y(t) can be very accu- 
rately approximated by y^ where 

Yg = l+-|p - /(l+iP)^-l 

and hence 3(t) and a(t) can be approximated by the following 
expressions . 

6(t) 2A(2+p-2y^) = (C.12) 

a(t) - 2A(y^^-y^(2-p)+l) = (C.13) 

Referring to Figure 5, we can see that when z is very large, 
such that the area under the curve 3(x) up to time x=1.0 
hour, is insignificant when compared to the area under the 
curve up to time x=z, then 3^z is a good approximation for 
z 

/ 3(x)dx. 

0 
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z 

2/3(x)dx 23z 

Similarly a(z)e ° can be approximated by a^e ® 

Thus an approximate solution for T(t) is 

a -23 t 

T(t) - Trf— [1 - e ® ] t > 10 hours 



However, 

T(t) 



for such large value of 




2(2+p-2y^) 



-23 t 

t, e ® ^ 0, 

t > 10 hours 



so we have 



58 



I 



I 



APPENDIX D 



Notations used in Simulation Program I 

a^ arrival time of next message at machine r 

departure time of next message at machine r, 

S =0 when machine r is not in use 
r 

denotes the machine to which machine r is transmitting 
N^ number of messages undergoing transmission at epoch 

CLOCK current clock time of simulation model 
e exponential variate with mean = 1 

X arrival rate of messages at each machine 

p departure rate of a message undergoing transmission 

N total number of machines 

n discrete uniform (0,N) variate 

TSTOP length of one simulation run 



59 




60 



I 



I 



I 






I 





61 



ono ooo ooooooooooooooooooooooonooooooooooooooooooooo 



SIMULATION 'LOSS' MODEL 
PROGRAM I 



SIMULATE A SYSTEM OF N. SENDOX MACHINES 
ASSUMING THAT ANY ARRIVAL THAT FAILS 
TO GET ON THE SYSTEM IS LOST 



A(J) 

S(J) 

LINK( J) 



NSYS 

N 

M 

TIME 

OELT 

P4 

P5 

PDF4 

PDFS 

ARRM 

DEPM 

TSTOP 

NRUN 

IX 

CLOCK 

NARR 

NON 



TIME OF NEXT ARRIVAL AT MACHINE J 

TIME OF NEXT DEPARTURE AT MACHINE J 

= 0 IF MACHINE IS NOT IN USE 

CONTAIN THE NUMBER OF THE -MACHINE THAT 

IS RECEIVING A MESSAGE OR TRANSMITTING 

TO MACHINE J 

NUMBER OF BUSY MACHINES 

TOTAL NUMBER OF MACHINES 

MAXIMUM NUMBER OF MESSAGES THAT CAN 

BE TRANSMITTED , = N/2 

RECORDS THE INSTANT OF ARRIVAL OR 

DEPARTURE 

DISCRETE TIME INTERVAL AT WHICH 

RELATIVE FREQUENCY ARE COLLECTED 

TIME AVERAGE PROBABILITY 

DISCRETE TIME PROBABILITY 

CUMULATIVE PROBABILITY FOR P4 

CUMULATIVE PROBABILITY FOR P5 

RECIPROCAL OF ARRIVAL PATE OF MESSAGE AT 

EACH MACHINE IN UNITS OF MINUTES 

RECIPROCAL OF DEPARTURE RATE CF 

TRANSMITTED MESSAGE IN UfJITS OF MINUTES 

LENGTH OF EACH SIMULATION RUN IN MINUTES 

NUMBER OF SIMULATION RUNS 

RANDOM NUMBER SEED 

CURRENT CLOCK TIME 

NUMBER OF MESSAGES THAT ARRIVE 

NUMBER OF MESSAGES THAT ARRIVE AND ARE 

TRANSMITTED 



DIMENSION A(200), S(200), LINK(200) 

2, NSYS(5000), TIME(5000) 

3, PKlOl), P2(101), P3(101), Y(lOl) 

4, DELT(IOOO) 

5, P4(101), P5(101) 

6t PDF4(101), PDFSdOl) 

CALL OVFLCW 

INITIALIZATION 

READ (5,1000) NRUN 
WRITE ( 6 , 1000) NRUN 
READ (5,1000) N 
WRITE (6,1000) N 
1000 FORMAT (15) 

READ (5,2000) ARRM, DEPM, TSTOP 
WRITE (6,2000) A RR M , OE PM, TS TCP 
2000 FORMAT (3F10.2) 

Ml=N/2+l 
DC 10 1=1, Mi 
P4( I )=0.0 
P5( I )=0 .0 
10 CONTINUE 
IX=52841 

START OF SIMULATION RUNS 
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DO 95 IC=1,NRUN 

IX=IX+I C*246 

CL0CK = 0 .0 

NARR=0 

NGN=0 

JAD = 1 

NSYS( JAD) =0 
TIME( JAD)=CL0CK 

AT THE START OF EACH SIMULATION RUN SET ALL S(K)=0 
AND SET A(K) EQUAL TO TIME OF FIRT ARRIVAL 



DC 20 I=1,N 
S( I )=0 

CALL EXPON (IX,E,1) 

Ad ) = E=^ARRM 
20 CONTINUE 

DO 25 I=L,M1 
Pl( I) =0 .0 
P2( I J =0 .0 
25 CONTINUE 
30 CONTINUE 

SUBROUTINE ARNEXT SEARCHES FOR MIN A(K),K=1,N 
SUBROUTINE ARNEXT SEARCHES FOR MIN S ( K ) , K= 1 ,N , S ( K ) »#0 

CALL ARNEXT (A,K,N) 

CALL DPNEXT (S,LtN,SM) 

IF (A(K) .GE. SM) GO TO 50 
IF (A(K) .GE. TSTOP) GO TO 60 

THE NEXT EVENT IS AN ARRIVAL AT MACHINE K 
TEST TO SEE WHETHER IT WILL BE TRANSMITTED 

CLOCK=A(K) 

CALL EXPON dX,E,l) 

A(K) = A( K)+E=^^ARRM 
NARR=NARR+1 

IF (S(K) .GE. CLOCK) GO TO 30 
40 CONTINUE 

GENERATE DESTINATION MACHINE NDES = 1 , N , N D ES «!K 

CALL RANDOM (IX,Z,1) 

NDES=Z=i=N 

NDES=NDES+1 

IF (NDES .EQ. K) GO TO 40 
IF (S(NOES) .GE. CLOCK) GO TO 30 
CALL EXPON (IX,E,1) 

E = E=^DEPM 
S(NDES) =CLOCK+£ 

L INKINDES ) =K 
S (K)=CLOCK+E 
L INK(K) =NDES 
JAD=JAD+1 
NCN=N0N+1 

NSYS( JAD) =NSYS( J AD-1) +1 
TIME( JAD) =CLOCK 
GO TO 30 
50 CONTINUE 

THE NEXT EVENT IS A DEPARTURE AT MACHINE L 
AND MACHINE JL (=LINK(D) 

IF (SM .GE. TSTOP) GO TO 60 
JL = LINK(L ) 

SUL) =0 
CLOCK=S (L ) 

S(L)=0 

JAL=JAD+1 

NSYS( JAD) =MSYS( J AD-1) -1 



63 



I 

I 

r 



I 

I 

I 

I 



I 







ooo oooo ooo ono 



TIME( JAD)=CLOCK 
GG TO 30 
60 CONTINUE 

TIME( JAD+1)=TST0P 

COMPUTE TIME AVERAGE PROBABILITIES 

00 70 I=ltJAD 

1 I=NSYS ( I )+l 

P1(II)=P1(II)+TIME(I+1)-TIME(I) 

70 CONTINUE 

COMPUTE DISCRETE TIME PROBABILITIES 
JJ=1 

NSTGP=TST0P+0. 1 
DO 80 !=1,NST0P 
DELT( I ; =FL0AT( I ) 

75 CONTINUE 
JJ=JJ+1 

IF (DELT(I) .GT. TIME(JJ)) GO to 75 

KK=NSYS ( JJ )+l 
P2(KK)=P2(KK)+1.0 
80 CCNTIMUE 

DO 90 1=1, Ml 
Y( I )=FL0AT( I-l ) 

P1(I)=P1( I )/TSTQP 
P2( n = P2( I )/TSTGP 
P4( I ) =P4( I )+Pl( I ) 

P5( I ) =P5( I )+P2( I ) 

90 CONTINUE 
IX=52841 

95 CONTINUE 
SP4=0.0 
SP5=0. 0 
VP4=0.0 
VP5=0.0 

COMPUTE STEADY STATE MEAN, VARIANCE AND 
CUMULATIVE PRGBABI LITI ES USING P4 L P5 

DO 99 1=1, Ml 
P4m=P4( I )/NRUN 
P5(I)=P5(I )/NRUN 
SP4=SP4+P4( I )*( I-l ) 

SP5=SP5+P5 (I )«( I-l ) 

VP4=VP4+P4( I )«(I-1)*(I-1) 

VP5=VP5 + P5 (I )»!=( I-l )*( I-i) 

IF (I .NE. 1) GO TO 96 
P0F4( I ) =P4( I ) 

PDF5( I )=P5(I ) 

GO TO 9 7 

96 CONTINUE 

PDF4(I )=PDF4(I-1)+P4(I ) 

PDF5( I )=PDF5( I-l )+P5( I ) 

97 CONTINUE 
11=1-1 

WRITE (6,8200) II, P3 ( I) , P4(I), P5(I), PDF4(I) 

1, PDF5(I) 

8200 FORMAT { I 10, 5F12 . 4 , / ) 

99 CONTINUE 

WRITE (6,8800) 

8800 FORMAT (////) 

VP4=VP4-SP4=:‘=:‘2 

VP5=VP5-SP5«*2 

PRINTOUT OF RESULTS 

WRITE (6 ,7000 ) NRUN,SP4,VP4 
7000 FORMAT (10X,'TI/^E AVERAGE PROBABILITY WITH', 15, 

1 ' SIMULATION RUNS ',/, lOX, • MEAN = ' , F 10 . 4 , ' VAR 1 ANC E 



64 



J 

I 



I 



( 



2 F10.4,//) 

WRITE(6,8000) NRUN,SP5,VP5 

8000 FORMAT (lOX, 'DISCRETE TIME PROBABILITY WITH', 15, 

1 ’ SIMULATION RUNS', /,10X, 'MEAN = ',F10.4, 

2 • VARIANCE = « , F 1 0.4 , /// / ) 

CALL PLOTP (Y,P3,M1,0) 

WRITE (6,8800) 

CALL PLOTP <Y,P4,M1,0) 

WRITE (6,8800) 

CALL PLOTP (Y,P5,M1,0) 

WRITE (6,8800) 

STOP 

END 

SUBROUTINE ARNEXT (A,K,M) 

DIMENSION A(IOO) 

K = 1 

AM=100000,0 
DO 50 1=1, N 

IF (A( I ) .GE. AM) GO TO 50 
AM=A( I ) 

K=I 

50 CONTINUE 
RETURN 
END 

SUBROUTINE DPNEXT (S,L,N,SM) 

DIMENSION S(IOO) 

L = 1 

SM=2000000.0 
DO 50 1=1, N 

IF (S(I) .LE. 0.00001) GO TO 50 
IF (S( I ) .GE. SM) GO TO 50 
SM=S( I ) 

L = I 

50 CONTINUE 
RETURN 
END 
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Notations used in Simulation Program II 



a arrival rate of messages at each machine 

d departure rate of a message undergoing transmission 

N total number of machines 

M = N/2 

rate of entering state i 
rate of leaving state i 
1ST initial state at time 0 

nvimber of messages undergoing transmission at epoch 
e exponential variate with mean = 1 

q uniform (0,1) variate 

TSTOP length of one simulation run 
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SIMULATION ’LOSS' MODEL 
PROGRAM 1 I 



TO SIMULATE THE TRANSITION OF STATES IN A 
BIRTH AND DEATH PROCESS 



NRUN 

M 

RARP 

RDEP 

RL( J) 
KM( J) 

IX 

TSTOP 
NSYS( J ) 

1ST 

PI, P4 

P2, P5 

PDF4 

PDFS 



NUMBER OF SIMULATION RUNS 
MAXIMUM NUMBER OF MESSAGES THAT 
CAN BE TRANSMITTED 

RATE OF ARRIVAL OF MESSAGES AT EACH 
MACHINE 

RATE OF DEPARTURE OF A MESSAGE 
UNDERGOING TRANSMISSION 
RATE OF ENTERING STATE J 
RATE OF LEAVING STATEJ 
RANDOM NUMBER SEED 

LENGTH OF EACH SIMULATION RUN IN MINUTES 
NUMBER OF MESSAGES UNDERGOING 
TRANSMISSION AT EPOCH T(J) 

INITIAL STARTING STATE 

STEADY STATE TIME AVERAGE PROBABILITIES 
FOR 1, 200 SIMULATION RUNS RESPECTIVELY 
STEADY STATE DISCRETE TIME PROBABILITIES 
FOR 1, 200 SIMULATION RUNS RESPECTIVELY 
CUMULATIVE STEADY STATE PROBABILITIES P4 
CUMULATIVE STEADY STATE PROBABILITIES P5 



DIMENSION RL(lOl), RM(lOl) 

1, NSYS(5000), T(5000I 

2, PKlOl), P2(101), P3(101), X(lOl) 

3, DELT(IOOO) 

4, PDF4(101), PDF5(101) 

5, P4(101), P5(101) 

CALL OVFLOW 

INITIAL IZATION 

READ (5,1000) NRUN 
WRITE(6,1000) NRUN 
READ (5,10U0) M 
WRITE(6,1000) M 
READ (5,1000) 1ST 
WRITE (6 , 1000 ) 1ST 
1000 FORMAT (15) 

READ (5,2000) RARR, RDEP, TSTOP 
WRITE(6 ,2000) RARR, RDEP, TSTOP 
2000 FORMAT (3F10.4J 
RAFR=1. 0/RARR 
RDEP=1 .0/RDEP 
M1=M+1 
M2=2«M 

RARR=RARPV (M2-1 ) 

P4( 1) =0 .0 
P5(l)=0.0 

COMPUTE RATES OF ENTERING AND 
LEAVING THE VARIOUS STATES 

DC 10 I =2, Ml 
I2=2«( 1-2) 

RL( I-l) =(M2-I2)=!'RARR*(M2-I2-1) 

RM( I) =( I-l )»i^kDEP 
P4( I ) =0 .0 
P5( I ) =0 .0 
10 CONTINUE 
I X=3564 1 
C 



69 



ooo ooo ooo ooo oooo oooo oooo ooo 



START OF SIMULATION RUNS 
REPEATED NRUN TIMES 

DO 95 NR=1,NRUN 
IX=IX+NR=^-513 
DO 15 1=1, Ml 
P1(I)=0.0 
P2( I ) =0 .0 
15 CONTINUE 

AT THE BEGINNING OF EACH RUN SET TIME T(0)=0 
AND NSYS(0)=IST 

IND = 1 

T(IND)=0.0 
MSYS( IN0)=IST 
20 CONTINUE 

CALL EXPON ( IX, E, 1) 

J=NSYS( IND)+1 
IF (J .EQ. 1) GO TO 30 
IF (J .EQ. Ml) GO TO 40 
RATE=RL(J)+RM( J) 

TdND+l) DENOTES THE INSTANT WHEN THE PROCESS 
CHANGES STATE 

T( IND + 1 ) = T( IND)+E/RATE 
IF (T(IND+1) .GE. TSTOP) GO TO 70 

TO DETERMINE WHETHER PROCESS GOES TO STATE 
J+1 OR J-1 

PB=RL( J )/RATE 

CALL RANDOM (IX,Y,1) 

IF (Y . LE. PB) GO TO 50 
GO TO 60 
30 CONTINUE 

PROCESS IS IN STATE 0 ALWAYS GOES TO STATE 1 

T( IND+1 )=T( IND)+E/RL( J) 

IF (T(IND+1) .GE. TSTOP) GO TO 70 
GO TO 50 
40 CONTINUE 

PROCESS IS IN STATE M ALWAYS GOES TO STATE M-1 

T(IN0+1 )=T(IND)+E/RM(J) 

IF (T(IND+1) .GE. TSTOP) GO TO 70 
GO TO 60 
50 CONTINUE 

PROCESS GOES TO STATE J+1 

K=NSYS( IMD) + 1 
IND=IND+1 
NSYS( IND)=K 

P1(K)=P1(K)+T( IND)-T( lND-1) 

GO TO 20 
60 CONTINUE 

PROCESS GOES TO STATE J-1 

K=NSYS( IND)-1 
1ND=IND+1 
NSYS( IND) =K 

PKK+2) =Pl(K+2) + T( IND )-T( IND-1 ) 

GC TO 20 
70 CONTINUE 

K = NSYS( IN0) + 1 
P1(K)=P1(K)+TST0P-T(IN0) 
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COMPUTE P2 THE DISCRETE TIME PROBABILITIES 

NST0P=TST0P+0. 1 
JK=1 

DO 80 I =1, NSTOP 
DELT( I J =I 
75 CONTINUE 
JK=JK+1 

•IF (DELT(I) .GE. T(JK)) GO TO 75 
JK=JK-i 
KK=NSYS ( JK)+i 
P2(KK)=P2( KKI+1 
80 CONTINUE 

DO 90 I =1, Ml 
X(I )=I-i 

P1(I)=P1( I )/TSTGP 
P2( I) =P2m/TST0P 
P4( I ) = P4( I ) + Pi( I ) 

P5(I)=P5( I ) + P2( I ) 

90 CONTINUE 
IX=35641 

95 CONTINUE 
SP4=0.0 
SP5=0. 0 
VP4=0.0 
VP5=0.C 

Calculate steady state mean, variance and 

CUMULATIVE PROBABILITIES 

DC 99 1=1, Ml 
P4( I )=P4( I )/NRUN 
P5(I)=P5( I )/NRUN 
SP4=SP4 + P4( I )*( I-I) 

SP5=SP5+P5(I )*(I-1) 

VP4=VP4+P4 ( I )*( I-l J-( I-l) 

VP5=VP5+P5 ( I )*( I-1)*(I-1) 

IF (I .NE. 1 ) GO TO 96 
PDF4( I ) =P4( I ) 

PDF5( I )=P5( I ) 

GO TO 97 

96 CONTINUE 

PDF4( I ) =PDF4( I-l )+P4( I ) 

PDF5( I } =PDF5( I-l )+P5( I ) 

97 CONTINUE 
I 1 = 1-1 

WRITE (6,8200) II, P3 ( I) , P4(I), P5(I), PDF4(I) 

1, PDF5(I) 

8200 FORMAT ( 1 10, 5F 12 . 4, / ) 

99 CONTINUE 

PRINTOUT OF RESULTS 

WRITE (6,8300) 

8800 FORMAT (////) 

VP4=VP4-SP4*«2 
VP5=VP5-SP5**2 
WRITE (6, 7000) NRUN,SP4,VP4 
7000 FORMAT (10X,‘TIME AVERAGE PROBABILITY WITH',!5, 

1 ' SIMULATION RUNS ',/, lOX, ' MEAN = ' , F 1 0 . 4 , ’ VAR I ANC E 

2 F10.4,////) 

ViR I TE( 6,8000) NRUN,SP5,VP5 

8000 FORMAT (lOX, 'DISCRETE TIME PROBABILITY WITH', 15, 

1 ' SIMULATION RUNS' ,/,10X, 'MEAN = ',F10.4, 

2 ' VARIANCE = ' , F 1 0.4 , //// ) 

CALL PLOTP (X,P3,M1,0) 

WRITE (6,8600) 

CALL PLOTP (X,P4,M1,0) 

WRITE (6,8800) 

CALL PLOTP (X,P5,M1,0) 

STOP 

END 
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THIS PROGRAM CALCULATES THE STEADY STATE ME AN , VAR 1 ANCE 
AND CUMULATIVE PROBABILITIES OF A BIRTH AND DEATH 
PROCESS 



M 

P 

CL 

CM 

PDF 

XMEAN 

XVAR 



DIMENSION X(201), P(201)i P0F(201) 

READ(5,1000) NCARO 
FORMAT (I5,2F10.4) 

DC 10 K=1,NCARD 

READ (5,1000) M, CL, CM 

HRITE(6,1000) M,CL,CH 

P(l)=1.0 

SUM =P( 1) 

M1=M+1 

M2=2*M 

CL=CL/(M2-1) 

DC 20 I =2, Ml 
J2 = 2=!=( 1-2 ) 

C1=(M2-J2)*CL*(M2-J2-1 ) 

C2=( I-l )>!=CM 

P( I) = P( I-1)''^C1/C2 

SLM=SUM+P( I) 

CONTINUE 
XMEAN=0 
XVAR=0. 0 
WRITE (6,2000) 

FORMAT (/////) 

DO 60 I =1,M1 
P( U = P( I ) /SUM 
X( I )=I-1 

XMEAN=XMEAN+P( I )«(!-!) 

XVAR=XVAR+P( I )*( I-l)*( I-l) 

IF (I .NE. 1) GO TO 30 
PDF(I )=P( I ) 

GO TO 40 
CONTINUE 

PCF(I)=PDF(I-1)+P(I) 

CONTINUE 

11=1-1 

WRITE (6,3000) II, P(I), POF(I) 

FORMAT ( I 10,2F15.4,/) 

CONTINUE 

XVAR=XVAR-XMEAN«*2 
WRITE (6,2000) 

WRITE (6,4000) XMEAN, XVAR 

FORMAT (10X,'MEAN = ' , F 10 .4 , lOX , • V AR I AN C E = * 
, F10.4,////) 

CALL PLOTP (X,P,M1,0) 

WRITE (6,2000) 

CONTINUE 

STOP 

END 



MAXIMUM NUMBER OF STATES 

VECTOR OF STEADY STATE PROBABILITIES 

ARRIVAL RATE OF MESSAGES AT EACH MACHINE 

DEPARTURE RATE OF A MESSAGE WHILE 

UNDERGOING TRANSMISSION 

CUMULATIVE STEADY STATE PROBABILITIES 

STEADY STATE MEAN 

STEADY STATE VARIANCE 
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Listing of GPSS Program 



This program simulates the arrivals and departures of 
messages at a system of 20 sendox machines. Functions EXP 
generates exponential variates with mean of one and UNI gen- 
erates discrete uniform variates between 1 and 20, The pro- 
gram consists of 20 closely similar segments. Each segment 
simulates the activities at one machine. At the end of each 
simulation run, this program gives the expected waiting 
times at each machine and steady state probabilities of 
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SIMULATION 'NO LOSS' MODEL 
GPSS PROGRAM 



C\J 

nO • • CVJOO 

• rvJrooo sO--4— 4 

o 

OOh-Ch 

ln^oc^c^ rOvOCT* 



^LU 

<2 


Chen 




—4 












*-hO 


000 C\J 


0 














r-H OC O 


LA • • r-4r^ 


AJ 


>- 












II <1-2 


• »-HrOr^ lA rH r— 4 


LL 


f— 
















+ 


»— 1 












< 




a> 














LU^00>- 


0 


»— 1 














:^oujcQ 


4*'00' lALALA 


UL 


0 












•-CNJ^ 


>^*OOQ>Ch (MIAC3 


4* 


< 




-0 








•“t — < H— 


• • • • • • • 


00 


LL 










’T> 


UJ.-4X^-H 




0 ^ 






cc 












-4LL 


1 — 




LU 






<0 


< <UJ 


lA 0> 


U.+ 


< 




GO 








hh2::s:q 


lAsOO^CNJ 


+ r- 






z 






—4 


o^ix: ^-huj 


ro • • • OvO 


C7^ —4 


to 




z> 








<o>“>5: 


• rHCSj'O »— 4 f— 4 


LLLL 


LU 




z 






to 


>UL to 




+ + 


0 










LU 


►-•rDCii- 




COnO 


< 




> 






»-H 


-J^CQ 


CO 


LL-4 


to 




H- 






1— 




UAO' 


+ U- 


to 




»— 1 








H-. UL20 


encoo'O^ cvjLAoo 


r- + 


LU 




-J 








1— UJQ<h-i 


• • • • • • • 


LLLA 


z 




*— f 






HH 


^1- to 




+ —4 






0 






0 


ujujoi^t/) 




sOLL 


0 




< 






< 


^oruj<>-H 


0000-4 


U- + 


z 




LL 






LL 


oocQ>5: 


Ajmeoro 


+ 4* 


»— 1 












CL to 5 : 0:00 


c\i • • • lA 


LA —4 


z 




Z 






X 


x»-«:dljj:^ 


• — 4CSJLA rO(J^'-4 


ULLL 


0 




0 






1- 


ujQ2:i— < 




+ + 


0 










0 


zee: 




'<f ro 


z 


Cl 


1- 






CQ 


UJ LU LU I— 


LA 


LL-4 0 


1— « 


X 


< -4 








l-l-X 


IA4-0 LALALA 


+ llajaj 




LU 


z > 






QC 


<<|— UJ 5 : 


(\jr*^o^G' — 4->^r— 


rO + 


LL 




►— « OC 






LU 


exTa: 1-0 




LLAJ-*.-)^ 


0 


z 


H- Zh- 






X 


UJUJUJOQ 


(NJ (\J 


+ rH AJQ. 




LL 


t/) ro »» 






\- 


ZZI-ZZ 


0 0 


AJLL>X 


-J 


•> 


»— 1 —4 






LU 


UJLUO'-<< 


esj »» 


LL + + UJ 


< 


0 


Q Zii^ 






X 


oozs:a: 


AJOOsjLAsOlA 


'+* —4 —4 


> 


N^ 


r— ILL <►* 






2 




5:-4 • • <}-0-4-4>2 




vt* 


•'—4 






LU 




a: •— 4og>t q:ov|oO'-4ajllulwul 


ca 


—4 


0-4 0. 




— 4AJ 


LU— J— 1 












Ql- 


LU 




LU CD 








< 




Z 


0 




to< 






LU 








LU 




-J 




Z OsJQN^ 


LULULU-JUJ 


LU 


LU 


u.< 


Z> 




O^-H 




0 »-4 ch O' 0 -4 


• _J t ICQh- 


1- 


K— 


0 LU 


0 




i-<o 




• • • **-H • • •— 4 QcAcO<< 


< 


< 


oozz 






>z 




1- 1- 


<<<^-H_| 


-J 


QC 


1— LUO 


CC 


LULU 


x< 


UJ > 


0 0 




z> 


LU 


f-H|— 


LU 


00 


0 LU 


Q_»-irsJ o: 


Z LA Z 


ce:a:a:<z 


z 


z 


>i-toto 


I- 


LULU 


LULUI— 


XZ-H UJ 


Z) r-nrocho 


<<<>HH 


f-H 


LU 


o:ol/>lu 


z 


Z>Z> 


Xo:< 


OJOtO to 


LL 0 > t *UL rOC^>>>U-LO 


L/> 


0 


l-Z<H 


LU 


CJfsJ 


oco 




0 *ojrA pHr^f— 4— 4 


















•UJ> 






—4 






—4 




a. »-H 


rsinc 






> 






-j 




X GOZLALALALA h-iUJ 






CC 






-j 




lu sOO'CMDOrOvOch— 4 cMtoi^n 






1- 






< 




0 • • • • • • 




* -«■ 






* •)? 



74 



CO 

o 

O 

O 

UJ 

CO 



a: 

UJ 





CO 




2 


















o 














oC 


















O 




00 














LL 




00 


















►— < 


Q- 




CL 








3 




5: 


X fvj 




X ro 






1— 1 






00 


LU > 


<NJ 


UJ > 


ro 










z: 




—1 




-J 




«J 






< 




-J 


22H 


-J 




< 


•-I 


> 


QC 


LL3 - 


< > 


0.3 •• 


< 


> 


»• 




OC 


K 




► a: 




•* 


q: 


•- 


CO 


m 




02X 


». m 


02X 


•* 


m 


s: 


m 


CO 


m 


stLL •« 


s: CO 


M-LL •« 


:s: 


00 


1— ii— 1 




f-H i-H 


O 


»*i— 1 


f— ♦ f— ♦ r—t 


•or-M 


f— <HH 


.-1 ^.-1 


Q-CO 


K 




< 


1— l»-I^CL 


m Nt CM CL CO OsJ Q- ro > CL fVj st 


.-1 I-H I-H CL m sO ro CL 00 rn CL vn > CL ro vO •-< 








00 














-J 




00 














»— * 




m 


m 




m 




UJ 


OC 


o 




5: 


H-UJ 


CJO 


HOI 


OC 


h- 


3m 


< 


mmm 




<»- m 


33 m mmm 


<H 01 


33 m 


mmm < 


z:u- 


LL 


KOCOcOK 


LL 




2:2Ll i^ococoh- 




220. 


HUC0C0H2 


CO 




mma:2:<<a: 


O 


hhCKIO I 


mm commo:2:<<cy:*-«Q:o mm 


'/> m moc 2 < < or: »-• 


m2: 


Q 


rsiM<<mm< 






ID 3 m u j 2 M < < m m < s: m • -H :3 3 m L u 2 


rM rvl< < m LU <I 2 


i-< 


m 


Q_>_j«ja 


Q 


nc: CO t/) u j m h- < H-i KH CL > m -j Q- o' 2 : CO CO m m h- h- <t 


•-•►-HCL>m_JCLa: 


<io: 


O 


mmmQmmm 


2 


mmcomiD3<<ciiujmmommmmmc;m33<<oL 


mmmQmujoim 




X 


C0C0Q<U1CCQ 


m 


H 0 < 1— 0 a <3 0 H - 1 / ) CO a <i q: cr. 0 h - 0 <x h- cj 0 0 H- 


cocoQ<ia:a.Qi— 



CM 




ro 


ro 


> 


-J 


>- 


-J 


on 




QC 


-J 


H 


< 


H 


< 



■Jt -K- -«•*■«■ 



75 



CL 






CL 






CL 






X ^ 






X lA 






X o 






UJ > 






UJ >- 


lA 




UJ > 


*o 






-1 




4^ < q: 


-J 










22H- 


— j 




22 h- 


U 




22h- 


— i 




LL3 •• 


< 


> 


LL3 •• 


< 


> 


ULO »* 


< 


> 




«• 


cC 


•*40 lA 




a: 


•*40 vO 


»» 


a: 




•< 


UJ 


02 it: 


«* 


lU 


02i«: 


»« 


UJ 


^Li. - 




00 


>tU- •• 


2 


to 


4-U. •* 


2 


00 


^ *<rH 


r— < HH f-H 




*J" •*«— 1 o 


f—i f— < r-H 


40i-H 


O '>J‘ •‘f — 4 » — i C\i 


r— 4 » — 1 


rH 1 — 1 40 C\J 


1 

1 

P 

7 

0 

4 
P 

5 
4 
P 

7 
V 
P 

4 

8 
1 
1 
1 
P 
9 
1 

5 
P 
S 
5 
P 
9 


V 
p 

5 
1 
1 
1 
1 
p 
1 
1 

6 
P 

5 

6 
P 
1 

V 
p 
6 
1 
1 








UJ 






UJ 




UJ 


UJ 


oC 




^-UJ 


oC 




»-UJ 


QC 


1 — 


K UJ 


OIDUJ 


UJUJIJU 


<H UJ 


IDOUJ 


UJUJUJ 


<K UJ 


IDIDUJ 


UJUJUJ < 


<22 


22U- h* 


OlOOOK 


2<22 


2211. 1- 


OtOtOH-2<22 


22U. 


HOOOkOH-2 


OCO UJIL) 


touJUJa: 


2<<o:»-«a:o ujuj 


ooujujq: 


2<<oc»-«a:e) UJUJ 


toujuj a: 2 < •* t oc 



OJ -H H O Z) UJ *J J 2: (Vl< < lU U 1< 21 Lii 1^. 3 D uu OJ 2: M rsj < < UJ UU < ^ OJ •-< 3 D : L) UJ ^ M rvl< < LU LU c ^ 

2 00 1/> Lu u j f- H- <i ^ a. > a. ft: 2: 00 oou./ uj 1 — < CL > _j a. a: 2 CO L u LU I— H a. > _j CL cr: 

UJ to UJ Z) ID < < D: U J UJ UJ O tU U J OJ UU UJ to UJ Z) 3 <5. < Q: UJ UJ UJ O UJ UJUJ UJ UJ OO U J 3 < < Q: UJ UJ UJ O UJ LJLJ LU UJ 



0<J. 1-0000}- 


toto Q < a: oc o h- o < K ooo o to to o <i oc or. Q h- o <i 


1— CJ50CJ 


-r 


< 


LA 


tA 


>o 


'O 


> 


— } 


> 




> 


-J 


oc 


—1 


oc 


_} 


a: 


-J 


1- 


< 


h“* 


< 


»- 


< 



76 



1 



I 



I 

. 5 * 



I 

I 






a. 






a. 






CL 






X f^ 






X 00 






X O' 






LU >- 


r- 




LU > 


00 




LU >- 


O' 






3 






3 






3 






3 




'z.'z.y- 


3 




z:^h- 


3 




LL3 •• 


< 


> 


U-3 


< 


> 


LL3 


< 


> 




•> 


OC 


•*4/^00 


•> 


CC 


•‘4^0' 


•> 


(X 




•> 


LU 




•> 


3 






3 


>tu_ •« 




t/) 


'tUL •« 


s: 


oo 


^UL •« 


5: 


KTi 


•»»— irO >4“ 


#— 4 ►— 1 


rO ^ ^ '4' 


•‘*—1 LTV \0 


«— IHH 


* — 4 LT\ -tA 


\0 >4“ ih^OO 




f— 1 4/^ 



00 

.-I —I Q_ r-^ h- Q_ 00 CL *-« > Q_ —I *-1 rH —I Q_ .-I —t ao O. CO 00 Q_ > Q- CO — I .-4 f-H Q_ .-I r-l Q_ (/ ) O' Q_ .-I > Q_ cr» —I 



LU LU 

Lu a: K-LU a: i — lu cc 

h- UJ 33UJ LULULU <1— LU 33UJ UJUJUJ <H LU 33UJ LULULU 

2:^u- h-oooooh-:z<f2^:^ z:^u- ^z:u_ ^-oa)oo^- 

OCO LULU 00LULUQCZ:<<CX:i-^QCO LULU OOLULUCrC ^<<QC*-hQCO UJLU t/)LULUOl <Q0 
LU i-< h- 3 3 lU LU 2: fvj M < < LU ai < s: LU H- 3 3 LU tU z: M M < < LU LU <1 5: OJ »-• H- Z) 3 LLI LU ^ M r>sj <r <l LU dJ < 

00 00 1 U u J f— K- < CL > _J Q_ or t/) t/) lU LU h- h- < KH KH CL > ^ Q_ ^ OO OO u J LU K K- <l — * hh CL > _J o. 

LU 00 LU 3 3 <l< QC LU LU LU O LU LU LU U J LU 00 LU 3 3 < <l Ct: LU LU LU O LU LU LU LU LU O') LU 3 3 <l< (X LU LU UJ O LU UJ LU 
0<I h-OCJOOH-00 00a<l0la:Ql-(3<K-U>CJOO|-i/)00O<CL'CXOHO<lKa><3OOt-0000Q*q:0r.Q:QH- 



> 3 

QC 3 

h- < 



00 

> 

CL 



00 



< 




77 



ERMINATE 



1 



I 

I 



I 



I 

I 



o 




r-4 






eg 






Q. »— < 




a. f-H 






Q. 






X > 


0 


X > 


*-H 




X > 


eg 




LU oc: 


f-H 


LU CC 






LU OC 


I— < 






_J 




-J 






— j 




2 : 2 ; V. 


-J 


2:2 


-J 




22 - 


-j 




U-3)0 


< 


> LUlDr-l 


< 


> 


UL30sJ 


< 


> 


r — 1 


•> 


QC r— I 


»» 


a: 


► r~H 


»> 


a 


0:2:^ 




LU 02^^ 




LU 


02 : 1 ^: 


»> 


LU 


>tLL 


2 : 


00 NfO. •* 




U1 


nMjL 




00 


•* r — 1 0^ C5 r— 1 f\| r-^ r— l»-H ^ 




g 'g’ ** 1 — ( cO'^ eg 1 


egi-iro4/)'r-i(\j >j’ 


Q. <>J rW Q. 00 CL > Q. (\J rH rW rH CL (\| CNJ r-l CL 00 CL (\l > Q- r-l (\J rH »~J r-l Q, (\J f\l CL 00 r-l a. f\l > CL »~j (\J r-< 






LU 






LU 




LU 


LU 


q: 


h-LU 


a: 




H-LU 


CC 


h- 


h- LU 


i:)3Lu 


LULULU LU 


i:)3Lu 


LULULU 


<h- LU 


ZDZDIU 


LULULU < 


< 2 : 2 : 


22U-. 




22LL. 


^-O0000^-2<t22 


22LL. 


H-O^/OLOl-2 


CjCO lulu 


L/)LUIUQ:2:<<QC»“^OCO lulu 


O0LULUQC2<<O:»-hqCO lulu 


l/) lU LU Qc: 2 < < a ►-< 



UJ D i:) UJ UJ 2 : rsj rvi < < ULI UJ< ^ UJ D 3 LU UJ rvj < < UJ ULI< UJ 1 ^. 3 Z) UJ UJ ^ fvl f < < Uj Uj < ^ 
^ oo uo ULi Ljj H* CL > _j u CL Li: 00 1/) Lu uj H* < HH »-i CL > _j _j cx 00 00 liJ uj H- H- < HH CL > _j ^ a. a: 

m 00 JJ II) 3 <x < q: uj m m Q Lu uj uj LU Lij c/) oj 2 ) 3 < < a: uj iJJ u j Q uj dj uj tu ULt oo uj II) 3 < <. (X uj ?xi uj o uj ULJ oj uj 
Lo<I^-c:?oe30^”^}ooQ<Ia.Q:a^-o<lH-c^a^oe)^—oo{/)a< Q:Q:o^—o<^^-c:>c?o(i)^—ooooQ<la:ol q>— 



o 

> 

a: 

h- 



o 



< 



> 

cc 



(M 



CNj 



>• -J 

-J 

h- < 



78 



I 



I 



I 

1 

I 

I 



cn 

CL 1 — 1 

X >- 


CO 




UJ CL 


r — 1 






_J 




zz: - 






u-3m 


< 


> 




•• 


CL 




•» 


LU 


4-LL •• 


X 


C/) 



•4’ •*f~4insOcri.— 1»— «rn.~<Lni^^rr) 
I-I I-I CL CNJ (\J I-I Q- 1/> ^ CL f\J > CL ^ 



4- 

O- r— 1 

X >- 


4- 




LU CL 








3 






3 




U-34* 


< 


> 


•*</> »— 1 


»> 


qL 


OX!iL 


»> 


UJ 


4*11. •• 


X 


00 



r-l I-I — HX OsJ Csi •"-! CL 00 ^ CL OsJ > CL I-I CsJ 



m 

CL r— 1 

X >- 


in 




LU CL 


r— I 






3 




zz: 


3 




U- 31 A 


< 


> 


f — 1 




CL 


oxx 


»> 


UJ 


stLL •« 


X 


00 



4' »*r— I r— I LO ^ r— llA O 

r-l CL C4 cn r-l CL 00 r-l X (\J > X r-l m 



UJ OJ 

uj a: LU c< j— U-) a: 

h- UJ Z>Z)UJ UJUJLU <f— UJ OOUJ UJLUUJ UJ 331X1 UJUJUJ 

2zu_ l-ooooo^-2<z^^ :^:?u- ^-ooooo^-:^<::^2: ^-ooooof- 

Qco ujuj ooLULijc^:z<<a:^-*a:o lulu uiuiuja:^<<.cjC^a:o ujlu ooujluq:^<<q(: 

LU K-^ h- c:) 3 LU LU rsi PnJ < ' i LU LU dil LU *-H J o UJ LU 2 : M rsj <X <t LU LU < 5: LU K 3 3 LU LU ^ r4 rsj < < UJ LU < 
71U) C/iLUtUI-l— <*'<»-iCL>_J-JCLQ::r'^00aJLUH-l— <»-'»-'0->3-JCLQ:2 :oOOOLULUI— I— <*-h^^CL>_J-JQl 
UJ o^; UJ 3 3 < < CL LU LU LU O LU LU LU U J UJ 00 UJ 3 3 < < CL LU LU I U O LU LU LU LU UJ c/ > LU 3 3 <l< tt: UJ LU » JJ o LU UJ L U 
o<i~aicju)cDj-ooooo<cLcLOH-o<i»~cjc?e)oi— ootj^o<cLQ:oKe;<u~o(jioOH*^i/)o<Q:cL:o. 




79 



TERMINATE 



GENERATE 1440,FN$EXP 
TRY16 ASSIGN 1»FN$UNI 

TEST NE P1,K16,TRY16 



vO 





r- 






CO 








Cl rH 






O. 3 








X >- 


r- 




X > 


00 






LU CL 


r-H 




LU cL 


r-4 








3 




40 >-H ' 


3 






22 •- 


3 




22 •- 


3 




> 


LL3r- 


< 


> 


LL3CO 


< 


> 


CL 


r— < 


•» 


CL 


r— < 


»> 


CL 


LU 


02: 




UJ 


02X 


»<• 


LU 


00 


nJ-U. 


2 


oo 


<tiL 


2 


00 



f\J vO 

rorOrHO-c/) 



%OrHr-<^r~lvOCNJ 4* ►rH CO 4* ^ rH xj" >t •'rH lA nO OQ r-* CQ lA ^ rH 00 sQ 

rH Q- ro > Cl rO rH a m ro rH Ql t/) r-H CL m > Q_ r-H CO rH r-( rH CL ro ro rH Ql C/) rH Q_ rO > CL .-H O'! r 



LU LU 

a: (— uj a: h-lu a: 

LULULU <H“ LU Z^Z^UJ LULULU <H UJ ID^LU LULUIU 

:^:^LL hztll h-oc/)ooK:?r<:z’ 2 : kooolo 

LULU c/)LutUA:2<<QiHHQ:o LULU ooLULua;2<<f cl:^-^q:o lulu ooluluoC2:<< 

3 LU LU :2 M M < < LU UJ X UJ *“• H - 3 3 LU LU 2 : M M < < LU LU < s: LU I — 3 70 UJ LU >I M isj < < L U LU 
LULUH-K-<-^^-«OL>3-Ja fi:2:t/>t;^miLUH-H-<--*^-4CL>_J_ICLQ::ZLOOOUI'>UHI— <--»^-«CL>33 
3 3 c< CL LU LU LU Q LU LU LU LU UJ oO LU 3 3 < < £X LU LU LU Q LU LU LU LU LU 00 LU 3 3 < < OC LU LU LU Q LU LU LU LU 
c:^c^OLD^^tA/o3Q<a:^xQ^-o<^-c300o^-o^c/)Q<cr:cLQKO<x^-c?<^oo^-~ooooQ<Q:^cah- 



^o 



< 



r- 

r— < 
> 
a: 
I-- 



r- 



< 



00 

rH 

> 

DL 

h- 



CO 



< 



80 



PART 

RMINATE 



Ch 






0 








CL rH 






CL <M 








X > 


o> 




X >- 


0 






LU £X 


r-^ 




LU CC 


<M 




rsj 


»— 1 J— 


-J 




♦ — ♦ 1— ^ 


-J 




UJ rH 


2 : 2 : 


-J 




2 : 2 : 


-J 




fsi 


LLZ>0' 


< 


> 


a. IDO 


< 


> 


►H rH 




»» 


q: 


*-40C\| 




oc 


00 


02X 




UJ 






LU 


^ 0 




5 : 


UO 


4-LL •- 


ST 


LO 


> *-0 


>X" *^rH CO 0 ^ O' r-H 4 /^ r-H 0^ CO 


'4' •'^H O' ^D CD rH »— 1 CD rH O' 4^^ rH CD CD CD 


•'rH »H0 


rH rH CL ro m --H cn/i r-^ X ro >« CL rH ro r-i r-H r-< X ro '4- c\i CL 00 c\i CL ro > CL c^j >x- — < 0 


rH OCLfNj 








LU 




LU 


LU 


LU 


UL 




I--IU 


DC 


KOJ 


LUH* 


K LU 


IDIDLU 


lululu 


LU 


IDIDLU 


LULULU <lh- 




<^2: 




^“OL000^- 


2:<2:2: 


2!2:ll 


KO0000H2:< 


2:<2: 



0^0 LULU i/)LULUQ:2<t<ta:^-^ci:o lulu LOLULuo:2:<<a:^-HL^e)-j*-'Lu»—H- 

UJ ID 3 LU lU r>J rvJ < <LU LU LU h- 3 3 LU LU M r>J.< < LU LU LU ID 51 LU 

2 1>0 00 LU LU H* IX > — > ^ CL q:: IT I/) oo ULT LL» t— J— < *-H CL > _J CL fx: OO 00 a: OQ < l>0 

LU O) LU Z) ID < <r OC LU IIJ LU a LU LU LU lU UJOO u J ID ID <1< OC LU LU LU Q LU LU ULI LU LU oO <l LU <. H- UJ 

<D< K cjf CP o (D H- Q <i oc cx. o o -X c/ o CD o H- i>o t/i a < u: o: Q exL t— K oo ct: 



0 


O' 


0 


0 


rH 


iH 


04 


04 


> 


-J 


D- 


-J 


QC 


-J 


ClC 


-J 




< 


1 - 


< 



81 



START 32000 

END 



APPENDIX E 



Calculation of X' 

In the previous sections, we had shown that a "no loss" 

model with input parameters N, A and p can be approximated 

by a "loss" model with input parameters N, X' and y, where 

X' is chosen such that the effective arrival rate X'(l-y^)^ 

for the "loss" model, is equal to X. 

When X.,(t) is in steady state, the mean number of busy 
M 

machines is Ny^, hence the probability of finding a busy 
machine is Ny^/N or y^. In the "loss" model, a message is 
transmitted only if it finds both the transmitting and re- 
ceiving machines free. The probability of this event is 
approximately assuming independence between avail- 

ability of machines and a large N. Therefore, the effective 
arrival rate, X for a "loss" model is X'(l-y')^ where 

y^ is given by equation (4.15) as 

y^ = l+-|p' - /(l+ip')^-l 

with p' = 2 ^ 

Therefore, we require that 

X' (/(l+^p* )2-l - ip')^ = X (E.l) 

After expanding the quadratic term in equation (E.l) and re- 
arranging terms, we obtain 

d+iP') - TTXr = /P' d+ip' ) (E.2) 

A p 
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Squaring both sides of equation (E.2) and recalling that 
p= — , we obtain 



p. = p + - - 2 



Making the substitution for p and p', we finally have 



X' = 



(1- i)" 

P 



Using this rate of arrivals in our "loss" model we are able 
to approximate the system size characteristics of the "no 
loss" model with arrival rate X. 
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